Evidence for X-ray synchrotron emission from simultaneous 
mid-IR to X-ray observations of a strong Sgr A* flare 
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ABSTRACT 

This paper reports measurements of Sgr A* made with NACO in L'-band 
(3.80 /zm), Ks-band (2.12 fim) and H-band (1.66 fim) and with VISIR in N-band 
(11.88 /xm) at the ESQ VLlQ, as well as with XMM-Newton at X-ray (2-10 keV) 
wavelengths. On 4 April, 2007, a very bright flare was observed from Sgr A* 
simultaneously at L'-band and X-ray wavelengths. No emission was detected 
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using VISIR. The resulting SED has a blue slope (/3 > for vLi, oc , consistent 
with vLy oc z/°'^) between 12 micron and 3.8 micron. 

For the first time our high quality data allow a detailed comparison of infrared 
and X-ray light curves with a resolution of a few minutes. The IR and X-ray flares 
are simultaneous to within 3 minutes. However the IR flare lasts significantly 
longer than the X-ray flare (both before and after the X-ray peak) and prominent 
substructures in the 3.8 micron light curve are clearly not seen in the X-ray data. 
From the shortest timescale variations in the L'-band lightcurve we find that the 
flaring region must be no more than 1.2 Rs in size. 

The high X-ray to infrared flux ratio, blue vL^^ slope MIR to L'-band, and 
the soft vLi, spectral index of the X-ray flare together place strong constraints on 
possible flare emission mechanisms. We find that it is quantitatively difficult to 
explain this bright X-ray flare with inverse Compton processes. A synchrotron 
emission scenario from an electron distribution with a cooling break is a more 
viable scenario. 

Subject headings: accretion, accretion disks — black hole physics — infrared: 
general — radiation mechanisms: general — Galaxy: center — X-rays: general 



1. Introduction 



The radio source coincident with the gravitati onal c enter of the Milky Way, named 
Sgr A*, was first discovere d by Balick & Brown in 11974 It had already been suggested 
( iLynden-Bell fc Reed Il97ll ) that the Milky Way may host a supermassive black hole at 
its center, and the newly discovered, unresolved source looked like it could well be the 
manifestation of such an object. That there really is a supermassive black hole of ~ 4 x 
IO^Mq, has now been proven beyond reasonable doubt through long-term monitoring and 
observation of the cluster of stars or biting within arcseconds of the black hole, most notably 
the star S2 (SO-2 in lGhez et al.ll2003l ) which has corapleted a complete 15- year orbit since the 
first monitoring observations in 1992 (jSchodel et al.ll2002l ; lGhez et al.ll2003l ; lEisenhauer et al. 
2005h . 



Sgr A* is thus a source of intense observational and theoretical interest, since it provides 
an avenue by which to study the physics of accretion in the presenc e of extreme gravitationa l 
fields. However, it is unusually dim for a supermassive black hole (IRieke &: Lebofskylll982l ). 
The spectral energy distribution of the radio source rises from radio towards submm wave- 
lengths, but no steady emission can be detected above roughly 10^^ Hz, implying that the 
spectral energy distribution (SED) turns abruptly around at this point (this feature has been 
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named the 'submm bump'). The overall luminosity is far below (by a huge factor of 
that expected for a black hole accreting at the Eddington rate. 



10^ 



It was only recently that the source was discovered at all on the high fr equency side of 



the s ubmm bump, where it was found to exhibit strong f lares in the X-ray (IBaganoff et al. 



2OOII ) and in the near- infrared (NIR) JCenzel et al.lboosh. A steady quiescent state in the 
X-rays at very low luminosities was also found ( iBaganoff et al.ll2003l ). The quiescent state 
has never been detected unambiguously in the NIR, nor has it ever been detected at mid- 
infrared (MIR) wavelengths for which only upper limits ca n be determined on either the 
quiescent state or possible flaring activity (see for example, ISchodel et al.l 120071 ). 



Subsequent to the first detections of Sgr A* flaring in the X-ray and NIR, a number of 
flares have been observed in both IR and X-ray wavelengths. Multiwavelength campaigns 
co-ordinating telescopes across the electromagnetic spectrum have worked towards obtaining 
simultaneous observations. 

Some general properties concerning the IR and X-ray flares that have emerged from 
those studies are: 



1. IR/NI R flares occur on average ~ 4 times per d ay (see e.g. Figure 18. lEckart et al. 
2006a ). or between 30-40% of the observing time ( Yusef-Zadeh et al. 20061 ). 



Strong X-ray flares occur on average ~ 1 per day (IBaganoffl l2003l ) . However, an en- 
hanced rate of X-ray flaring can be observed within a time interval of rou ghly half a day 



e.g. a bright flare followed by three flares of more moderate amplitude iPorquet et al. 



2008|). 



3. 



Every X-ray flare appears to be associated with a NIR flare, how ever not every NIR 
flare is associated with an X-ray flare (e.g., iHornstein et al.l 120071 ). 



4. X-ray and NIR flares occur simulta neously, with no significant delay (lEckart et al. 



2004J . l2006d : lYusef-Zadeh et al.ll2006[ ) 



5. Substructural variations with characteristic timescales of 15-25 minutes are seen in IR 



flares on a regular basis (IGenzel et al.l l2003l : iMeyer et al.l l2006bl : lEckart et al.l |2006c 



Trippe et al.ll2007h . 



6. Significant drops in flux are sometimes seen during X-ray flares (IBaganoff et al.ll2003 



Porquet et al.ll2003r ). 



7. Polarimetric investigat ions of the 



significantly polarized (lEckart et al 



ares in the NIR have shown that the source is 



in the tail end of the flare (iTrippe et al.l 120071 : iMeyer et al.ll2006bl ) 



2006b[) and that the polarizat ion angle can swing 
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8. At high fluxes t he flare has a constant blue spectral index i n ul]^ of /9 = 0.4 between 
3.8 and 1.6 fim (IHornstein et al.l 120071 : iGillessen et al.l 120061). For low fluxes it appears 



that the flare shows red uL^, spectral indices ( iGhez et al 



2005: Eisenhauer et al 



Krabbe et al.ll2006l ) with a poss ible trend of spectral index with flux (IGillessen et al. 



20061 ) although this is disputed ( Hornstein et al. 2007 ). 



2005 



The two brightest X-ray flares (jPorquet et al.l |2003| . |2008| ) have been observed to have 
well constrained soft uL^, spectral i ndex values P = 0.2 ± 0.3 and /3 = 0.3 ± 0.3 cal- 
culated at a 90% confidence range (jPorquet et al.ll2008l ). While several fainter flares 
were observed, only a small number of photo n index values has b een reported; the 



latter exhibiting harder spectral indices (e. g.. iBaganoff et al. 



200ll ). A re- analysis of 



XMM-Newton archival flares performed by iPorquet et al.l (120081 ) with a homogeneous 
data analysis shows that at low X-ray flux the spectral index is in fact not well con- 
strained and a soft index as found for the two brightest flares cannot be excluded. 
Similarly, Mascetti et al. 2008 (submitted) analyses a co-added spectrum of all Chan- 
dra flares to date and reaches the same conclusion (i.e. that soft uL,^ indices are not 
excluded). Therefore, higher S/N spectra for individual weak/moderate X-ray flares 
are still required to establish whether all flares have similar spectral shape or not. 

10. The X-ray flares appear unambiguously to be 'events', i.e. s hort, large amphtude 



outbursts followed by what looks like a perfectly flat baseline ( iBaganoff et al.l 12003 



Porquet et al.ll2008l ). In the infrared, it is less clear whether this picture applies or 



whether the IR 'flares' are simply peaks within an underlying sea of variability with 
the characteristics of red noise. Similarly, it is debated whether the substructural 
features seen in IR flares correspond to a characteristic frequency of the system (a 
quasi periodic oscillation (QPO)), or whether it is caused by statistical fluctuations in 
a smooth, red noise power spectrum (T iMeyer et al.ll2008l ). 



Note that there are also many important results from observations at longer wavelengths, 
but since we are directly concerned with the IR and X-ray data we have obtained, we have not 
gone into them in this paper. The full results of our April 2007 multiwavelength campaign 
including the observations at radio and submm wavelengths will be presented in Yusef-Zadeh 
et al. (2009, in prep.). 

The quiescent state of Sgr A* can be successf ully described by either a radiatively 
inefficient accretion flow ( or RIAF: see for examp le lYuan et al.ll2003l ). or as arising from 
the base of a compact jet (IFalcke fc Markoflll2000l ). Each of these models can describe the 



^Here and elsewhere in this paper we use (3 to denote the i^L^, spectral index, defined as z^X^ cx 
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observed properties of the quiescent state with similar magnetic field strengths {B ^ 30 
G) and electron energies (7 = E/mc^ ^ 10). The addition of Bremstrahlung emission 
from within the Bondi accretion radius (-RBondi ~ 1") explains the X-ray quiescent emission 
jQuataertlbooi ). 



The origin of the fiare emission within either of these basic pictures is much less certain. 
The high degree of linear polarization of the flares at IR wavelengths points to a synchrotron 
origin, but the emission mechanism responsible for the X-ray flares is not known. In the 



Eckart et al. 


(2004 


) and 


, . ^^^^ 

Eckart et al. 


f 2006c) 



explained the simultaneity and the observed fluxes of infrared and X-ray flares through 
the synchrotron self Compton (SSC) emission of a compact source component emitting pri- 
marily at mm/submm wavelengths, with the emissi on at IR waveleng t hs po ssibly due to 



a combination of synchrotron and SSC emission. In lYusef-Zadeh et al.l (120061 ). the IR and 
X-ray observations were interpreted within a picture where the X-ray emission was due to 
inverse Compton scattering of subm m and IR p h otons i nvolv ing populati o ns of both submm- 



emitting and IR-emitting electrons. iLiu et al.l (l2006al ) and lYuan et al.l (l2003l ) also present 



models involving i nverse Compton proc esses. Synchrotron mod els for the X-ray flare have 



been sugge sted bv iMarkoff et al.l (1200 ll ) and I Yuan et al.l (120031 ) and IR/X-ray synchrotron 
models by Yuan et al. (j2004 ). Synchrotron models in general have been criticized due to 
the fact that the high energy electrons needed to generate X-ray synchrotron emission have 
very short cooling timescales (much shorter than the typical X-ray flare duration), requiring 
continuous injection in order to replenish the high energy population. However, this may 
not be such a disadvantage, and continuous injection is in fact a natural and reasonable 
expectation for the kinds of processes responsible for particle acceleration such as magnetic 
reconnection, turbulence and shocks. 

Going beyond the emission process behind the flare, there are models which attempt to 
simultaneously describe the detailed properties at one wavelength , such as the hot spot model 



(IBroderick &: Loeb 



2005 



models (iTagger &: Melia 



Meyer et al. 



2006a 



20061 : iFalanga et al.ll2008r ) 



Trippe et al.l 120071 : ?), or accretion instability 



In this paper we present our multiwavelength observations and focus on constraining the 
emission mechanisms responsible for the simultaneous IR and X-ray flares we have observed. 
Although a full analysis of the detailed time-resolved SED evolution is beyond the scope of 
this paper, our high quality, full coverage, lightcurves in both L'-band and X-ray wavelengths 
offer the for the flrst time the opportunity to undertake detailed modelling of the time 
evolution of the flare SED, which may shed further light on the emission mechanisms and 
physical conditions/processes giving rise to a flare event. 
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^Throughout this paper we adopt a Galac tic Center distanc e of 8 kpc (lEisenhauer et al. 

20031 ). and a black hole mass of 4 x 10^ Mq (T IGhez et al.ll2008l ) for which the Schwarzschild 
radius is Rs = 1-2 x 10^^ cm. For the solar luminosity we used the value Lq = 3.8 x 10^^ 



erg s 



2. Observations 

In this section we present IR/NIR (3.8, 2.1 and 1.6 /im), MIR (11.88 yum) and X-ray 
(2-10 keV) observations of Sgr A* taken in April 2007. In particular we focus on April 4, 
2007, on which date a very bright flare was observed in both L'-band (3.8 /im) and X-ray 
simultaneous to the MIR observations. 



2.1. IR/NIR Observations 



The IR/NIR observations were taken at the VLT in Chile as part of a multiwave- 
l ength campaign (LP 179.B-0261) in A pril 2007, using the NAOS-CONICA instrument 
( iLenzen et al.ll2003l : iRousset et al.ll2003l ) in imaging and polarimetric modes. We observed 
between 5:00 and 11:00 UT on April 1 to April 6 obtaining data in L' (3.8 /im), Ks (2.1 /xm) 
and H (1.6 /zm) wavelength bands. 

We subjected the raw data to a sky subtraction computed from jittered object images in 
the L'-band case, and from dedicated observations of a patch of sky devoid of stars ~ 700" W 
and 400" N of the GC for the Ks and H band observations. This was followed by fiat-fielding 
and a correction for dead/hot pixels. 

Once we had reduced the set of images, the raw fiux at the position of Sgr A*§| in each 
image was determined via two independent methods: (i) aperture photometry, where the 
fiux was computed as the sum of all pixels within a small aperture centered on Sgr A*, from 
which the sum of pixels (normalized by area) within a larger annular region surrounding Sgr 
A* was subtracted t o remove backgroun d contamination; and (ii) PSF photometry, where 
we used StarFinder (jPiolaiti et al.l l2000l ) to automatically identify and extract PSFs from 
the reduced images, thereby obtaining source fluxes. 

Finally, we calibrated the raw fiux with the fluxes of nearby stars of known and stable 
brightness, and converted it to a physical fiux. For the extinction correction we used the 



■^confused in these observations with the star S17. 



- 7- 



values Al = 1.8, Ak = 2.8 and Ah = 4.3 flGenzel et al.ll2003h . 



For those observations taken in polarimetric mode, we added the fluxes obtained in 
ordinary and extr aordinary images to obtain an integrated flux (for further details see 



Trippe et al.ll2007h . 



The resulting lightcurves for the combined source Sgr A* + S17 are presented in Figure 
dJ Several weak flares are seen (labelled 1-4 and 6-7). On the night of 4 April and under 
good conditions (seeing ~ 0.55 — 0.9 and Strehl ratios ~ 0.45 — 0.65), a very strong flare 
was seen in L'-band at the position of Sgr A*, beginning just before 06:00 UT, April 4, and 
lasting for roughly 2 hours. We present the lightcurve of this flare in detail in Figure [2j 

Since Sgr A* was confused with S17 on April 4 2007, S17 also contributes flux to the 
lightcurve shown in Figure [21 In addition, the quiescent state of Sgr A*, if it exists, is not 
well known so it is possible that the qui escent state and po ssibly other L'-band sources such 
as a small dust cloud close to Sgr A* (IClenet et al.ll2005l ) also contribute to the minimum 
flux of Sgr A* on the night of April 4. It is also not clear from the lightcurve whether Sgr A* 
ever reaches a level of non-activity during our L'-band observations. If we take the mean of 
the group of points where the lowest flux for the night was recorded between t ~ 295 to 300 
minutes, this leads to an upper estimate for the contribution of any nonvariable emission of 
4.3 mJy. 

The combination of the data quality and the strength of the flare activity make the April 
4 flare presented here the best specimen of our entire dataset recording the L'-band activity 
of Sgr A* and spanning 2003 to 2007. The lightcurve shows very significant substructure 
on a timesc ale of ~ 20 rninutes. This kind of sub s tructure has been seen in previous Ks- 
band flares (iGenzel et al.ll2003l : lEckart et al.ll2006bl : iTrippe et al.l 120071 ). but is seen here for 
the first time in L'-band. This strengthens the case that the presence of substructure is a 
common feature of IR/NIR flares. 



2.1.1. Limits on L' -hand flaring activity ~ 7-11 hrs UT. 

On April 4 2007 we also obtained data in other wavelength bands (Ks-band polarimetry, 
H-band imaging and Ks-band imaging) for another ~ 3.5 hours after the flare observed in 
L'-band (Figure [T]). There was no obvious variability in these wavelength bands. We find 
that our measured Ks-band flux for Sgr A* + S17 of Sk = 5.1 ± 0.4 mJy is consistent 
with the measurement of Do et al (2008) for S17 alone (5.5 mJy, dereddened according to 
rriK = 2.8 as assumed in this paper; no error was given). It is also consistent with our own 
past measurements of S17's magnitude {rriK = 4.8 mJy), although our error on this value is 
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Fig. 1.— Lightcurve for Sgr A* + S17 for April 1-6 2007 IR/NIR observations with NACO 
at the VLT. Observations were taken at L' (black), Ks (red), H (blue) as well as in Ks- 
band using polarimetric mode (green). Several flare events are seen over the six nights of 
observations, labelled 1-7. Some nights show a more continuous level of variability (April 3), 
while on other nights there are long periods with no obvious variable emission. The L'-band 
flare from April 4 is the most significant event seen. A flare of equivalent strength in Ks 
band (given a colour of (3 = 0.4 in uL^) would reach Sk ~ 20 mJy. On April 5, another 
L'-band flare is seen under less favorable conditions which reaches Sl' ~ 20 mJy in L'-band. 







300 320 340 360 380 400 
time [minutes after 00:00 UT, April 4] 



420 



Fig. 2. — Extinction-corrected L'-band flux of Sgr A* on 4 April 2007, determined via PSF 
photometry (black data points) and via aperture photometry (grey data points). Above the 
lightcurve of Sgr A* is shown in red (PSF photometry) and light red (aperture photometry) 
the flux of the nearby star S2 as a flux comparison, clearly demonstrating that the substruc- 
ture is intrinsic to the Sgr A* source. S2 was confused at the time of observation with S13, 
and in the flgure their combined flux is shifted upwards by 25 mJy. Sgr A* was confused 
with S17 and the flux shown also possibly has a contribution from a dust cloud as well as the 
unknown quiescent state of Sgr A*. The data are binned to a bin- width of 44.3s. The PSF 
photometry method only lists fluxes of sources detected with 3cr signiflcance, which explains 
why the lightcurve derived by PSF photometry begins only at around ~ 300 minutes. We use 
the mean of these points with t < 300 minutes as an upper estimate of the background level 
(S17 -|- confused sources -|- quiescent state) upon which the flare emission is superimposed. 
If the minimum of the lightcurve is used instead, we obtain a background estimate of 2.4 
mJy. 
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very large and on the order of 2 mJy. The H-K color of the combined Sgr A* + S17 source 
is P = 2.6 ± 0.7, consistent with that of a pure stellar source {(3 = 3.0). From this and 
the lack of significant variability in the emission we would conclude that in H and Ks band 
the emission is dominated by SI 7. We note that in deconvolved images of this dataset we 
do see an elongation of the source in Ks band, and resolve two distinct sources in H band. 
Of these two H-band sources we can not be sure whether the source coincident with Sgr 
A* is stellar (e.g. from faint unresolved S-stars surrounding the black hole) or whether it 
might be quiescent/flaring emission. Due to the close proximity of the sources (only 3 pixel 
separation) it was not possible to determine the Ks-band fluxes of each source accurately. 

In L'-band the lowest measured luminosity lies above the extrapolation of the (Sgr A* 
+ S17) H-K measured color and we can conclude that within the timespan of our L'-band 
observations we did not reach the flux level of S17, which would be expected at ~ 1 mJy. 
The remaining flux we see may be due to a further contamination of the L'-band flux by a 
confused source (e.g. a small dust cloud known to be an L'-band source nearly coincident 
with Sgr A*, ~ 4.7 mJy Clenet et al 2005. although this seems unlikely given the high flux), 
or it might also be due to the fact that the flare activity never ceased within our L'-band 
observation time interval. Since flares are redder than the stars, we can not rule out that 
some low-level flaring continued to occur after ~ 430 minutes while we observed in Ks- and 
H-bands. 

To make some estimation of the L'-band flux during the time interval within which we 
observed in Ks and H bands we must extrapolate from our Ks-band measurements, which 
introduces large uncertainties. We can reasonably assume that any flaring emission was 
below Sk ~ 2.1 mJy during the Ks and H-band observations (our lower limit on the flux 
of S17 is 3 mJy; note this is also consistent with extrapolating the flux of the deconvolved, 
separated source in the H-band images to Ks-band with a slope oi P — 3). Although the 
flare color at low flux levels is not well established, if we take a flare color of vL^ — 0.4 
we can estimate that the combined source of Sgr A* + S17 should not have been at a flux 
level higher than 4.7 mJy in L'-band during this time. However, if the flare were redder at 
low flux levels or if some flux is contributed from the dust cloud near Sgr A* then the limit 
on the flux level could be higher. We obtain an upper limit of 9 mJy if we add the lowest 
flux detected during the L'-band observations (4.3 mJy), using it as an upper limit on the 
magnitude of any nonvariable contamination. 
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2.2. X-ray Observations 



On the 4th of April 2007, the VLT observations described in Section 12.11 overlapped 
with those of XMM-Newton. The observations and data reduction of the X-ray flare are 
published by Porquet et al (2008), and will not be repeated here. The X-ray lightcurve is 
presented in Figure [31 where it is compared with the L'-band lightcurve. 

The X-ray flare was very bright. The 2-10 keV spectrum of the flare showed a soft 
spectrum: a power law fit correcting the underlying model for dust scattering and absorption 
gives a power law slope of F = 2.3 ± 0.3 (error bars given at the 90% confidence level), 
equivalent to a v Ly spectral index of j3 = — 0.3±0.3. The bright flare observed on April 4th 
(labelled #2 in iPorquet et al.ll2008l ) is the second brightest flare obs erved so far from Sgr 
A* with an amplitude of about 100 compared to the quiescent state. iPorguet et al.l (120081) 



show that this flare and the brightest X-ray flare observed in October 2002 (IPorquet et al. 



20031 ) have similar light curve shape, duration, and spectral characteristics (photon index). 



2.3. Mid-Infrared Observations 



VISIR, the VLT Imager and Spectrome ter for the mid-infrared, mounted on th e ESO/VLT 
telescope Melipal (UT3) at Paranal, Chile (iLagage et al.ll2004l : iPantin et al.ll2005l ). observed 
the Galactic Center from 2007-04-04 05:29:00 to 2007-04-04 10:34:00 (UT). We collected the 
data with the imaging PAH2_2 filter on, at 11.88 ± 0.37 /im in the atmospheric window "N". 
The Small Field mode (SF) was employed, resulting in a field of view of 256 x 256 pixels 
(19.2 arcsec^), each pixel corresponding to 0.075 arcsec^. 

We performed the calibration of the PAH_2 filter on 2007-04-04 05:16:24 fUT) with a 



109.9 s observation of the standard star HD 102461 (9.237 Jy in the PAH2_2 filter: ICohen et al. 
1999[ ). 



The basic 'chopping and nodding' technique was applied to acquire the data, which 
were then reduced with the standard VISIR pipelin^: this involved flatfielding, bad pixel 
correction and combination of a stack of chopped and nodded frames to produce a final set 
of 79 consecutive images. 

We determined the position of Sgr A* using the precise positions of the SiO maser 
sources IRS 7, IRS 9 and IRS lOEE, distributed about Sgr A* jReid et al.ll2007h . No point 
source at the position of Sgr A* is detected in either the individual images or the collapsed 



^see http : //www. eso . org/ instrument s/visir/ 
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Fig. 3. — Comparison of the L'-band and 2-10 keV lightcurves (labelled IR and X, respec- 
tively). The top two panels show the two flares over the period of the L'-band observations. 
In the lower two panels, a larger time interval is shown. We also show data taken in other 
wavelength bands (Ks-band polarimetry, H-band imaging, and Ks-band imaging) subsequent 
to the L'-band flare indicating that the flare activity ceased in both wavelength bands af- 
ter roughly 4:00 UT. In the same night of observations three more X-ray flares were seen 
(Porquet et al 2008); the first of these started at UT 11:32 (692 minutes), i.e. roughly an 
hour after the last of the NIR observations. The dashed lines indicate our estimates of the 
background levels (i.e. emission that is not flaring emission) at each wavelength. In the case 
of the IR lightcurve this may be an overestimate. 
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image of the entire night. We also performed a Lucy- Richardson deconvolution with HD 
102461 as PSF with again no source detection. The flux from a box of 0.375 arcsec^ centered 
on the position of Sgr A* is constant with an average value of 123 ± 6 mJy. This flux may be 



Eckart et al. 


2004; 


^ ^^^^^^ 
Schodel et al. 


2007) 



To determine an upper limit of the brightness of the simultaneously observed flare, we 
subtracted an average image of the quiescent phase (12 images from 07:17 to 08:00) from an 
average image of the flaring phase (12 images from 05:30 to 06:13@). We included a slight shift 
in the relative positions (less than half a pixel) of the images, degraded the mean images with 
Moffat functions to mimic the slight differences of atmospheric conditions between them, and 
finally destriped the subtracted image. Over a region of ~ 3 arcsec^ centered on Sgr A*, 
north of the minispiral, the subtracted image displays a relatively flat background. 

To quantitatively estimate our detection limit, we proceeded by simulating an artificial 
flare in the data. We included a weak point source of a given flux (with VISIR's PSF) at 
the position of Sgr A* in the substracted image. We increased the point source's flux until it 
was detected at a significance of 3a and took this value as an upper limit on the flare's mean 
flux. We thus estimate that Sgr A* could not have been brighter than ~ 12 mJy at 11.88 /im 
(3(7, not dereddened). Note that this value is compatible with VISIR's empirical sensitivity 
at this wavelength: 7 mJy/lOcr/1 hr (median value for different atmospheric conditions). 

The value of the extinction correction in the MIR depends critically on the strength and 
shape of the silicate absorption feature at ~ 9/im. The values in the literature are published 
as ratios relative to Ak or Ay, so we use the value Ak = 2.8 mag {Ay = 25) mag to ensure 
consistency across our multiwavelength observations. The closest extinc tion measurement 
to our observation wavelength, A = 11.88/im, was m ade bv iLuta (119991) fo r a wavelength 



of ^ 12.4/im. We consider three theoretical models ( Chiar fc Tielensr2006 : Draine fc Lee 



1984J : iRoche fc AitkenI Il984l ) for the shape of the silicate profile in the region, to allow us 
to extrapolate the value measured at 12.4/im to 11.88/im. These models each use different 
sources as template profiles but are all very similar in slope arou nd 12 fim and result in very 
similar values of extinction when normalized to the iLutzl (119991 ) value at 12.4 fim. The thus 
determined extinction value is An gSfim = 1-7 ± 0.2 mag. With this value, the dereddened 
3a upper limit on emission from Sgr A* during the flare is F^^^ ~ 57 mJy. 



^Note that VISIR observations started 5 min after the beginning of the X-ray flare, which was from 05:25 
to 06:13 
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3. Results 

3.1. Simultaneity of infrared and X-ray flare 

We have obtained complete, fully sampled lightcurves in L'-band and X-ray. The X-ray 
lightcurve showed no other events during either a ~ 15 hour period before and ~ 5 hours 
after the bright X-ray flare. In the IR/NIR, in the remaining observations of about 3| hours 
following the L'-band flare, there were no other flares or obvious variable emission above 
a level of Sk ~ 2.4 mJy (implying < 12% the equivalent Ks-band flux of the peak of the 
L'-band flare, using j3 = 0.4). We translated this to a more conservative limit (see Section 
I2.1.ip in L'-band of Sl' < 9 mJy (< 32% the peak L'-band flare flux). At least on this 
occasion, both infrared and X-ray emission can thus be best characterized as isolated 'flare' 
events. 

From a correlation analysis, the L'-band and X-ray flares are found to be simultaneous 
to within ~ 3 minutes. In particular we do not see any signiflcant delay or asymmetry in 
the longer wavelength emission relative to the peak of the X-ray flare, thus excluding that 
adiabatic expansion of an initially optically thick blob plays a role in the infrared and X-rays. 

3.2. General Lightcurve Shape 

This mult i wavelength observation allows us to make the most detailed lightcurve com- 
parison so far, of simultaneous IR and X-ray flares from Sgr A*. From the comparison of 
the two lightcurves shown in Figure [3] it appears that the L'-band flare begins flrst. The 
L'-band lightcurve rises before any significant X-ray emission is seen, and L'-band emission 
remains after the X-ray emission has subsided. It appears that the two events have different 
durations from one another, and that the infrared event lasts longer overall than the X-ray 
event. 

Taking the uncertainties on the background levels into account, we measure FWHM 
durations for each (background-subtracted) lightcurve of FWHM/zj = 66±8 and FWHMx = 
28 ± 0.5 minutes. Thus we find that the FWHM of the L-band fiare is ~ 2 times that of the 
X-ray fiare. 
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3.3. Substructure 

One very striking feature of the L'-band lightcurve is the substructural variations seen 
on a timescale of ~ 20 minutes. The variations in flux are large: up to ~ 30% the peak flux. 
There are no apparent features that would correspond to these in the simultaneous X-ray 
lightcurve. 

Although the X-ray lightcurve has lower SNR than the L'-band lightcurve, the error 
bars are only on the order of ~ 10% the peak flux, and thus cannot hide substructures as 
large as those in the L'-band lightcurve. From this we conclude that the lack of substructure 
in the X-ray lightcurve is not due to lower SNR, and that this property is in fact intrinsic 
to the simultaneous IR/X-ray lightcurves. 



3.4. Shortest time-scale variations 

In the L'-band lightcurve, in particular at t ~ 350 minutes but also near to t ~ 370 
and 395 minutes, very rapid changes in flux (factors 120% to 170%, significance > 3a) are 
observed within a very short timescale. At < 47 seconds. 

Such short term variations place a limit on the size of the flaring source, or at least 
the size of the part of the source providing the sudden change in flux (which is a significant 
fraction, ~ 30%, of the total flux). Since such variations cannot propagate within the source 
faster than the speed of light c, the source size Rp is immediately constrained to be 



Rf < cAt = 1.2Rs. 



A caveat to this constraint is that we have not considered various relativistic factors; 
given the small size obtained and the fact that we think the flare might occur at very small 
radii in the accretion flow, relativistic effects might be an important factor in influencing 
the time variability of the source. Relativistic beaming near the event horizon is a possible 
source of magnifying the amplitude of variations due to an underl ying spatial structure in the 
infrared emission. Such beaming effects could be considerable (T IBroderick fc Loebll2005l ). 



3.5. Power Spectra 



Whether or not the substructures seen in the L'-band lightcurve are indicative of a 
QPO or are merely spurious peaks in a red noise process is a matter of current debate 
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(e.g. T IMeyer et al.ll2008l ). Since the putative QPOs inevitably turn out to be too weak to 
stand a significance test from a single observation niglit's wortli of data, we turn our focus 
to longer timescales. H-band and Ks-band polarimetric data that were taken following the 
L'-band measurements show no evidence of any variable emission, and this clearly holds an 
implication for the variability behaviour of the source on longer timescales. 

Figure H] shows the periodogram of the L'-band data compared with that of the X-ray 
lightcurve. We use the ordin ary periodogram with the RMS-squared normalization (see. 



for example lUttley et al.l |2002| ) which allows us to compare lightcurves taken with different 
instruments (and here at different wavelengths). For a consistent comparison between the IR 
and X-ray variability we took the mean from the same time interval, i.e. the maximum time 
overlapping time interval of the IR and X-ray observations. For timescales < 130 minutes 
(frequencies > 0.008 min~^) we show the power spectrum of the L'-band data only. We 
use our limits on the variable emission in Ks and H bands to constrain the periodogram at 
lower frequencies. Some uncertainty in the normalization of the IR power spectrum comes 
about through our uncertainty in the mean value given the extrapolation from Ks-band to 
L'-band. There is an apparent peak at ~ 20 minute timescales. Whether it is a real QPO 
or the spurious peak of a red noise spectrum, it is noteworthy that the putative QPO peak 
of the L'-band data has no corresponding feature in the X-ray power spectrum. This is 
consistent with our observation that the substructures of the IR lightcurve are not present 
in the X-ray lightcurve, which is comparatively smooth. 

At low frequency we also see the difference in widths of our lightcurves; the power 
spectrum of the X-ray lightcurve resembles a Gaussian at low frequencies which is as expected 
for the power spectrum of a single Gaussian-like flare event. The power spectrum of the IR 
lightcurve resembles a narrower Gaussian, again expected from the fact that the IR lightcurve 
was of longer duration than the X-ray lightcurve. We note then that the clear flattening of 
the power spectrum towards low frequencies again suggests that the IR flares are discrete 
events. 



3.6. Spectral Energy Distribution 

The 'flare state' SED for the observations of Sgr A* on April 4, as determined by our 
multiwavelength observations, is shown in Figure O 

For the vLy value at L'-band we computed the mean of the extinction corrected, back- 
ground subtracted L'-band flux. We chose to take the mean value rather than the peak 
value since the MIR limit and X-ray spectra were both determined as averages over the flare 
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Fig. 4. — Power spectra of the L'-band (thick black) and X-ray (blue squares) lightcurves. 
Our constraint on the periodogram at low frequencies from Ks and H band data is shown 
as the gray region (with the interval mean shown as black dashed line). The RMS-squared 
normalization was used, where the mean of the IR and X-ray lightcurves was taken from the 
same time interval. A peak is seen around the 20 min timescale (~ 0.04 min^^), while there 
is no corresponding peak in the X-ray power spectrum. The two dashed lines indicate the 
corresponding FWHM frequency for the Fourier transform of a Gaussian given the FWHM 
durations of the IR and X-ray flares. 
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V (Hz) 

Fig. 5. — The Spectral Energy Distribution of Sgr A*: in black (filled triangles) ar e radio 
to submm measurements of the quiescent state (IMarkoff et al.ll200ll : IZhao et al.ll2003l ). Note 
that these measurements are time averaged measurements and the errorbars include variable 
emission of up to 50%. As open black circles wi t h arro ws are shown 30 /xm, 24.5/im and 
8.6 jU-ni upper limits taken from iMelia fc Falckd (120011 ). the upper limit at 8.6 /zm from 



Schodel et al.l (120071 ) and the limit on t he quiescent s tate a t 2 /im from iHornstein et al. 



(120021 ). The quiescent state values from iGenzel et al.l (120031 ) are shown as the open black 
circles with errorbars, and the flare values from the same paper as gray flUed circles. The 
X-ray quiescent state is shown as the black bow-tie ([Baganoff et al.l 120031 ) . The dashed line 
shows a model for the quiescent state ( Yuan et al. 20031 ). Our new measurements for the SED 
of a flaring state of Sgr A* are shown in red (fllled squares): (i) the MIR (11.88 /im) upper 
limit is shown as the downwards-pointing arrow. The MIR upper limit is determined over 
an interval 05:30 to 06:13 (see discussion in Section [2731) . (ii) The L'-band measurement is 
shown as the red square with errorbars. This corresponds to extinction corrected, background 
subtracted mean value of the L'-band observations, 19.1 ± 3.6 mJy (in this case we used the 
minimum of the lightcurve, 2.4 mJy as the background estimate). Since the MIR observations 
did not start until 5:30 (approximately half an hour after the onset of the NIR flare), the 
mean was computed over the MIR time interval rather than the entire L'-band flaring interval. 
Also shown next to the L'-band data point is how the L'-band measurement would continue 
into Ks and H band with a slope of 0.4, characteristic of the L'-H slope of a numb er of 



bright' observed flares in the literature (IHornstein et al.l 120071: iGillessen et al.l 1200 



]^ and 



also consistent with the slope of the peak flare values of iGenzel et al.l (120031 ). (iii) Two 
possible X-ray spectra are shown, neither of which is model-ir idependent. The red points 
which slope upwards to the left indicate the power law flt, of iPorquet et al. J2008I ). while 
the grey points show the blackbody flt of the same paper. The X-ray spectrum was scaled 
by a small factor since it incorporated data from an extra flve minutes before the MIR 
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interval. The error in the L'-band data point is computed as the standard deviation of the 
hght curve. 

For the X-ray data, it is difficult to show an intrinsic, dust and absorption-corrected 
X-ray spectrum without the assumption of a model. This is because the inversion of a raw 
counts/ channel X-ray spec trum is gener ally non- unique and unstable to small changes in the 
counts/channel spectrum (lArnaudlll996l ). To determine the best fitting spectrum, a model is 
calculated and 'folded', or convolved with the instrumental response after which the folded 
model spectrum is compared to the observed counts in each channel. Once one has found a 
best fit, the process can be reversed for the best fitting model and one obtains an intrinsic, 
but model dependent, spectrum. 

Because of this, we can not show a single, model-independent X-ray measurement on 
the SEP for the April 4 flare. Instead, we show two possible X-ray spectra as obtained 
by iPorquet et al.l (120081 ) (i) assuming a power law shape, with F = 2.3 ± 0.4 and Nh = 



12.8; ^;? X 10^^ cm-^ (at the 90% confidence level, using the statistic; see Appendix 



B in iPorquet et al.ll2008l ). and (ii) assuming a blackbody model, with parameters Nh = 
X 10^2 cm-2 and kT = 1.51°;} keV. The blackbody fit had the lowest A^^ of the models 
investigated in that paper. Both models show a soft spectral index above u ~ 10^^ Hz. We 
scaled the X-ray data, which was determined over the full X-ray flaring interval, by a factor 
0.95 to account for the fact that the interval over which the MIR upper limit was calculated 
was shorter by 5 minutes (the scaling factor was determined as the ratio in fluxes between 
these two intervals). 

In both L'-band and X-rays, there is a substantial increase in flux above the quiescent 
level. The absence of any detectable emission at 11.88 /im implies that the flare emission 
spectrum must ri se from 11.88 /iin to 3. 8 fim. This appears consistent with a (u L,,) spectral 
index of /5 = 0.4 jHornstein et allboOTl : Icillessen et aPbood : ICenzel et al]l2003h . 



The rise in uLi^ from the MIR towards NIR wavelengths suggests that the population of 
electrons producing the L'-band flare must have a different distribution of electron energies 
to those in the submm bump. This might be an power law tail of transiently accelerated 
electrons, for example, or a small group of electrons heated to a high temperature. What 
this observation shows is that a NIR flare cannot be due to a small change in overall prop- 
erties of the submm bump (such as, for example, a global increase in magnetic field which 
temporarily increases the emitted synchrotron emission of the quiescent state). The flare 
event must involve only a small fraction of the quiescent state electrons, either in some kind 
of acceleration process that acts globally but inefficiently within the accretion flow, or via a 
more efficient but very localized acceleration process. As mentioned in Section 13. 4[ the short 
time scale variability of the source also points towards a localized event. 
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4. Modelling the flare state SED of Sgr A* 

We studied the flare state SED of Sgr A* under four simple SED models. These models 
explore different emission scenarios for the IR/X-ray flare and correspond to scenarios where 
the IR and X-ray flares are due to synchrotron and inverse Compton emission mechanisms. 
For the fitting of the S ED models to the data we use the X-ray spectral fitting program 
XSPEC JArnaudI EqQgI i To incorporate our infrared data points we add an extra data 



channel with the L'-band extinction corrected fiux; the instrumental response for this data 
channel is an identity matrix. The four XSPEC models we used are: 

1. icmodel: IR emission is synchrotron emission by transiently heated/accelerated elec- 
trons; X-ray emission is due to the inverse Compton scattering of submm photons from 
the 'quiescent' population of electrons by the population of electrons producing the IR 
emission. 

The seed photon spectrum comes from the quiescent population of electr ons and its 



spect rum and total luminosity is fixed (we use the model spectrum of lYuan et al. 



20031 ). i?Q, the size of the region containing the quiescent state (submm-emitting 
electrons) electrons is a free parameter and controls the photon density of submm 
photons available for inverse Compton scattering. 

The IR synchrotron emission is modelled based on a thermal distribution of electrons. 
There are three parameters that pertain to the transiently heated population of elec- 
trons producing IR synchrotron emission: B, the magnetic field, 6e, the dimensionless 
electron temperature {6e denotes the typical energy 7 of the electron distribution; it 
is equal to kTe/mc^, where is the temperature of the thermal electron distribution), 
and N, the total number of IR synchrotron emitting electrons. 

2. sscmodel: IR emission is synchrotron emission by transiently heated/accelerated elec- 
trons; X-ray emission to IR/NIR photons of the transiently heated/accelerated (fiare) 
electron population that are inverse Compton scattered by the same population (i.e. 
SSC). 

The IR synchrotron emission is again modelled based on a thermal distribution of 
electrons. This IR synchrotron emission can be again computed from the parameters: 
B, magnetic field, 6e, the dimensionless electron temperature (the typical 7 of the 
electrons, see icmodel), and N, the total number of IR synchrotron emitting electrons. 

In this model it is the parameter Rp, the size of the region containing the fiaring 
(IR-emitting) electrons that controls the photon density of the seed photon spectrum. 
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3. powerlaw: IR emission is synchrotron emission from a power law energy distribution 
of accelerated electrons. The parameter of interest in this model is the particle index 
of the power law electron distribution, p, i.e. A^(7) oc 7"^. 

4. powerlawcool: IR emission is again synchrotron emission from an electron distribu- 
tion with continuous injection of power law electrons and the addition of synchrotron 
cooling. Parameters of this model are p, the particle index of the injected electron 
spectrum, and the magnetic field B, which determines the cooling time of electrons 
and thus the energy /frequency at which the cooling break occurs. 

The two models icmodel and sscmodel were developed especially for this work, while 
powerlawcool was a simple adaptation of the existing XSPEC model powerlaw to incorpo- 
rate the cooling break. Specific details of the icmodel, sscmodel and powerlawcool models 
are listed in the Appendix. 

For each model we also take into account the effect of photoelectric absorption and 
dust scatte ring on the X-ray spe ctrum via the XSPEC routines scatter and wabs (for more 



details, see iPorquet et al.ll2008l ). These effects were not applied to the NIR data. For the 
dust scattering (scatter), we fix Ay = 25 to match the dust extinction corrections used for 
the L'-band and MIR data. For photoelectric absorption (wabs) we allow the parameter Nh 
to be determined. 

As an extra constraint, we add an extra data point at H-band (1.65 fim) which cor- 
responds to a constraint on the uL^ slope from L'-H band of (3 = 0.4±0.2. We find that 
adding this constraint generally results also in models which do not violate the MIR 3a upper 
limit. We ran models also without this extra constraint and very similar best fit values were 
obtained. 

Table [T] lists the parameters of the fitted models. Figure [H] shows the SED corresponding 
to the best fitting case for each model. In the next subsections we go through each model 
in detail. 



4.1. icmodel: Flare caused by inverse comptonized submm bump photons 

A best fit model for the case of submm photons scattered by IR-emitting electrons is 
shown in Figure O The model is a satisfactory fit to the data. The typical electron energies 
involved appear reasonable (7 ~ 140). However, the magnetic field strength of 210 G is high 
compared to that expected for the inner regions of the accretion flow (10 — 30 G), and the 
parameter Rq has a best fit value of 0.046 Rs, which as we will argue is an unreasonably small 



Table 1. Models: Fit Parameters 



Parameter 


Icmodel 


Fit to Mean Fluxes & NIR spectral index 
sscmodel pouerlav 


powerlawcool 


Nh [x1022 cm-2] 


11.7 


(9.9, 14.3) 


11.5 


(9.7, 13.7) 


11.5 (10.6,12.7) 


12.4 


(11.0, 12.1) 


B [Gauss] 


210 


(30, 2900) 


6000 


(2200, 7900) 


< 0.1, or > 60 


6.1 


(0.1,60) 


Be [kTe/nieC^] 


140 


(50,210) 


11 


(9, 16) 








Ne [xlO*o electrons] 


4.7 


(0.2, 130) 


1.5 


(0.7,4.2) 








Rq [Rs] 


0.046 


(0.001,0.27) 












Rf [Rs\ 


> 0.02 




0.0013 


(0.0009, 0.0020) 








p 










2.88 (2.82, 2.94) 


2.4 


(2.1,3.1) 


/ d.o.f. 


70.1 /74 




69.9 /74 




72.4 /77 


70.4/76 




reduced x'^ 


0.95 




0.94 




0.94 


0.93 




Violates 3<t MIR upper limit? 


No 




No 




Yes 


No 





Note. — Summary of best fit parameters for different scenarios: synchrotron + submm IC (icmodel), synchrotron + NIR SSC (sscmodel), 
simple power law (powerlaw), power law with cooling breaJf (powerlawcool). In each case a NIR vLi, slope of /? = 0.4 ± 0.2 was enforced 
in order to add enough constraint to the parameters. We found that models which violated the NIR slope by > 2cr usually violated the 3it 
MIR limit also. Listed is also whether the model violates the MIR limit. Next to each value we provide the 90% confidence interval for each 
parameter. 
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Fig. 6. — The best fits of tlie four XSPEC models for the April 4 flare IR and X-ray data. 
The X-ray data points are the unfolded spectrum for the given model (blue solid line); only 
the PN unfolded spectra are shown, with the data points binned for plotting purposes, (i) 
icmodel: The solid blue line shows the best fit IC model for the fit constraints of the April 
4 flare, which satisfies the MIR limit and the NIR spectral index. The dashed blue line 
shows the best fit model holding the magnetic field fixed at B=30 G. This model violates 
the NIR spectral index and comes close to violating the MIR upper limit. It does however 
allow a larger size for the quiescent region, Rq = 0.27 Rs- This is however still far from 
the size/photon densities expected from size measurements of Sgr A*, (ii) sscmodel: Best 
fit SSC model. In this case the magnetic field and density are extremely high. The source 
becomes self-absorbed in the NIR, and the spectrum shows strong curvature from L' to H- 
band. (iii) powerlaw: Best fit power law model. This model violates both MIR limit and 
NIR spectral index, (iv) powerlawcool: A more feasible synchrotron model with a cooling 
break. This model corresponds to the steady state solution for a system with a constant 
injection of power law electrons where the energy loss of the electrons due to synchrotron 
emission is taken into account. 
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size to contain the quiescent state of Sgr A* (stated another way, it corresponds to a much 
higher photon density of submm photons within the flaring region than can be expected). 

To understand whether the small value of Rq is really 'too small', we must be sure of 
how the parameter Rq should be interpreted. We note that -Rq is in fact constrained through 
the quiescent photon energy density required to produce X-ray emission of the amplitude 
that we see in the X-ray flare. Since photon density is defined through Uph = L/cA, with A 
the surface area of the region emitting the luminosity L, we see that Rq can be interpreted 
as a constraint on the surface area of the (quiescent) emission region. 

Thus the true quiescent region's geometry must have a surface area equivalent to the 
surface area of a sphere of radius Rq in order to reproduce the required photon density. For 
example, a torus situated at the last stable orbit, -Rlso = S-R^ would have an equivalent 
surface area to our best fit Rq for a ring thickness of 2 x 10~^ Rs'- extremely thin. 



New observations of the size of Sgr A* at 1.3 mm, (IDoeleman et al.ll2008l ) approaching 
the peak of the submm bump, suggest that the quiescent emission region may not be centred 
on the black hole. The reason for this is the fact that the measured size for Sgr A* at 1.3mm 
is smaller than the minimum apparent (gravitationally lensed) size allowed for an object very 
near a black hole. If the quiescent emission region is indeed offset from the position of the 
black hole then we no longer require an extremely thin ring containing most of the quiescent 
region electrons. However, this still does not solve the size issue: at the 90% confidence 
level, the largest value of -Rq that is compatible with the data, -Rq = 0.27 R5 (0.54 in 
diameter) is still far below the measured FWHM size at 1.3 mm of ~ 3.7 R^. 

Another issue stems from the fact that our model does not take into account a stratified 



region (i.e. the p roperty that the observed size changes with wavelength, iBower et al.ll2004j : 
Shen et al.ll2005l ) in the calculation of the photon density from -Lsubmm and Rq. Within the 
true (non-homogeneous) quiescent source the local density of quiescent photons will change 
with position. It may be more realistic, rather than to input the photon density via the 
variable parameter -Rq, to input the known photon density of a known model for Sgr A* 
which satisfies all the observations including the size m easurements . As a demonstration of 



how we can implement this we again take the model of I Yuan et al.l (l2003l . see Figure E]) for 



which we have obtained tables of the quantity n(i/, -R) at different radii -R. Thi s model has 



been shown to predict sizes at 3.5 and 7 mm consistent with those observed (lYuan et al. 



20061). 



Implementing this model, the free parameter replacing -Rq is the radial position of the 
flare from the the central black hole, Vpos- This position will determine the local photon 
energy density that is to be inverse Compton scattered by the flare electron population. 
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When we incorporate the model's photon density into our inverse Compton code however, 
we find it difficult to find a reasonable fit for any Vpos- 

One can understand why this is through Figure [TJ which shows the energy density 
spectrum for (i) our simple 'one-size' models, for the best fit Rq as well as a ~ Sa upper 
limit (i.e. corresponding to the lowest energy density spectrum allowed to produce acceptable 
X-ray IC scattering) compared with (ii) the energy density spectra for various inner radii 
of the Yuan et al model. As we can see, the photon density in the model is just too low 
over the entire frequency range to reproduce a bright, soft. X-ray fiare via inverse Compton 
scattering. 

Finally, although we were fitting an IC model to the fiare, we must not forget about 
the fact that the electrons producing the NIR and IC emission must also be producing SSC 
emission. In fact, assuming the IC scenario to be the cause of the X-ray flare, we can put 
a lower limit on the size of the flare emission region Rp, by requiring the absence of an 
SSC contribution. To estimate this, we slowly varied Rp from its maximum value (-Rq) 
and observed at which point SSC emission began to overwhelm the IC X-ray emissioii^. We 
found that the flare region must be more than 0.02 Rs in size. This would imply a density 
of ng < 8 X lO^cm"^, which could be compatible with the kinds of densities (n^; ft i lO^cm"^) 



expected for the inner regions of an accretion flow near Sgr A* (lYuan et al.ll2003l ) 



4.2. sscmodel: Flare caused by inverse Comptonized NIR flare photons (SSC 

case) 

The critical frequency for synchrotron emission {i>c, the frequency at which a synchrotron- 
emitting electron emits most of its energy, and thus the energy at which the uLiy spectrum 
turns over, see Section 16. ip is linear in B and quadratic in 7. Thus to obtain an SSC peak 
below 10^® Hz the electrons must have very low energies of 7 ~ 10-15. If the electrons were 
to have such low energies, it follows that a very large B is needed to produce a synchrotron 
peak above NIR frequencies so that an increasing lyL^y spectral index is observed in the NIR. 

Thus it makes sense that our best flt SSC model, shown in Figure [6l corresponds to low 
electron energies and high magnetic fleld strength. The magnetic fleld strengths required are 
enormous, a factor of > 200 greater than the typical magnetic flelds of the quiescent state. 
At the same time, a dramatic decrease in the ratio of 6e/B will have the effect of suppressing 



^For the electron energies of our IC models, the SSC emission peak always occurred higher than X-ray 
frequencies, so the SSC spectrum had a hard vL^ spectral index and is not a valid solution. 
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Fig. 7. — Comparison of the photon density at each frequency in (i) the best fit IC model 
with Rq = 0.046.^ .^ ? (thic k red sohd hne) and (ii) different radii between 1.5 Rs and 10 Rs 
in the lYuan et al.l (120031 ) model (thin blue solid lines). The dashed red line shows the the 
photon density spectrui n for Rq = 0.27, a t the 90% confidence level for the parameter Rq. 
The photon densities of I Yuan et al.l (120031 ) are in general 2 to 3 orders of magnitude too low, 
which shows why there is no well fitting IC mode l for the mu l tiwave length observations of 
the April 4 flare, given the photon densities of the lYuan et al.l ( 12003! ) model. 
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the SSC emission. Accordingly the density required from the small size, rig ~ 9 x lO^^cm"^, is 
equivalent to a density enhancement on the order of 10^ above densities typical for the inner 
regions of the accretion flow and is also unrealistic. The high densities required even have the 
result that the synchrotron spectrum at IR/NIR wavelengths becomes self-absorbed. The 
self- absorption results in a very steep spectrum at IR wavelengths, and it shows significant 
curvature. Overall, due to the extreme physical conditions required to create the observed 
X-ray emission via SSC, we rule it out as the emission process behind the April 4 X-ray flare. 

4.3. powerlaw: X-ray and IR flare from single power law synchrotron emission 

We have found that neither the IC nor SSC scenarios are entirely satisfactory as explana- 
tions for the simultaneous observations of April 4. Therefore, it is worthwhile to investigate 
other possible scenarios for the production of the X-ray flare. One such possibility is that 
both IR and X-ray flares are synchrotron emission. 

A power law energy distribution of electrons 

N{E)dE oc E-'PdE 

in the presence of a magnetic field, will create a synchrotron emission spectrum 

pL, oc z/(=^-p)/^ 

While the powerlaw model gives a reasonable fit to the L'-band and X-ray data, it 
violates the MIR limit and gives a soft vLy spectral index in the NIR. This model is also 
unrealistic since for reasonable magnetic field strengths S ?a 10 — 30 G we expect electrons 
within the energy range of our power law distribution to have very short cooling timescales. 
Either very low {B < O.IG) magnetic fields are needed or very high {B > 60G, together 
with a very flat spectrum of injected electrons p ~ 1.9) to prevent a cooling break occurring 
between IR and X-ray wavelengths (i.e., this motivates our next model powerlawcool). We 
can thus definitively rule out this scenario. 

4.4. powerlawcool: X-ray and IR flare from power law synchrotron emission 

with cooling break 



It is well known for synchrotro n emission source s to exhibit various breaks in their 
spectra due to cooling processes (e.g. iPacholczykI Il970l ) . One of the lowest frequency breaks 
likely to occur is due to synchrotron losses. The electrons responsible for the emission above 
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this cooling break lose energy due to synchrotron cooling faster than they can typically 
escape (which they do on roughly the dynamical timescale). If the source of acceleration 
in the plasma occurs continuously (i.e. there is a continuous injection of electrons from the 
heating/acceleration process), a steady state solution exists where the spectrum follows the 
usual synchrotron spectral index of (3 = (3 — p)/2 (with p the particle index) at energies 
below a characteristic energy, the 'cooling break', while above this ene rgy the spectral in dex 



flattens to /? = (2 — p)/2, corresponding to a particle index of p + 1 (lYuan et al.ll2003l : ?). 
The position of the cooling break corresponds to the electron energy (or, in the emitted 
spectrum, the frequency) at which the cooling time is equal to the escape time. The escape 
time is uncertain but for our simple estimate we will use the dynamical time: 

B \ ( V \ -1/2 



Td,n = y^^5min (2) 

with R ^ 3.5 Rs the radial position within the accretion flow. The cooling break therefore 
occurs at a frequency of 

Vcooi = 2.56 X 10^^ Hz (3) 

This model has more freedom than the power law model, and the data provide less 
constraint on physical parameters than in the IC or SSC (because in those cases the IR and 
X-ray flares arise from different emission mechanisms). We can determine what magnetic 
fleld strengths are necessary for such a model since the magnetic fleld B directly influences 
the position of the cooling break. 

Such a cooling break model flts the data well. The magnetic fleld strengths we flnd for 
this case are of the order of the magnitude of those expected for the inner regions around 
Sgr A*. 



5. Flare Evolution: Lightcurve Shape and Substructure 

The SED modelling we presented in the last section only examined the mean proper- 
ties of the flare emission. Our observations hold a great deal more valuable information in 
the time-dependent properties of the lightcurves. It is interesting to examine what different 
emission scenarios imply for the evolution of simultaneous flare in the IR and X-ray bands. 
There are two outstanding features of the simultaneous lightcurves that need to be under- 
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stood: the broadness of the NIR hghtcurve in comparison to the X-ray hghtcurve, and the 
substructures seen in the NIR hghtcurve but not in the X-ray hghtcurve. 

The synchrotron and IC luminosities depend on 5, 6e (think 7)0, and the size of the 
quiescent region Rq as 

Lsynch (X NelB'' (4) 
Lie cx NOIRq^ (5) 
The SSC luminosity depends not on Rq but on the size of the flaring region Rp as 

Lssc oc NH%B^R-/. (6) 

The (changing) parameters governing the overall flare evolution are probably N{t) and 
the electron temperature 9E{t). For the case of submm photon IC, these are the only pa- 
rameters that can affect both the optically thin synchrotron luminosity (not dependent on 
flare or quiescent region size) and the IC scattered luminosity (not dependent on magnetic 
field). 

For the submm IC case, both the synchrotron and inverse Compton luminosities depend 
on the same powers of N and 9e- This means that if only N or 6e were to change throughout 
a flare, the X-ray hghtcurve should follow the same functional form as the synchrotron 
hghtcurve. That is, both lightcurves should have the same width, or duration. We can take 
this either 

• as further evidence against the submm photon IC scenario, or 

• to imply that if the X-ray flare is to be explained by the submm IC scenario, then the 
observed hghtcurve widths can only be produced if some parameter aside from N and 
9e also varies throughout the flare. 

If we consider the second possibility, then the magnetic field, B, is perhaps the most 
obvious choice for the varying parameter. The dependence of the lightcurves on B means 
that for a change in B to create a broader NIR hghtcurve, the magnetic field must decrease 
during the flare. The process must reverse itself towards the end of the flare: the magnetic 
field must increase again towards initial values. Such behavior could, for instance, occur 
if stored magnetic energy in a small region were released to accelerate electrons, as in a 



^We use 9e and not 7 in our arguments because 9e represents a characteristic energy of the entire 
population of electrons, while 7 more properly denotes the energy of each electron in the population. 
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reconnection event. Remember however the high magnetic fields required {B ~ 200 G) for 
the submm IC picture with icmodel in Section |H which makes it difficult to decrease the 
magnetic field during a flare unless we accept even higher values for the magnetic field before 
and after the flare. 

The only other possibility to explain the lightcurve durations in the submm IC scenario 
is that the photon density increases during the flare (i.e. effectively through the parameter 
Rq). Although it is not realistic for the overall photon density of quiescent state photons to 
change much with time, the photon density experienced by the ffare electrons could increase 
as it moves inwards within the accretion ffow towards higher submm photon density. Again, 
to explain the second half of the ffare, for this possibility the position of the ffare within the 
accretion ffow must move outwards towards lower photon density to increase the IR emission 
relative to X-ray emission as the ffare declines. 

As far as the substructure is concerned we can see from the same relations that a 
variation in the magnetic ffeld affects the synchrotron luminosity but not the IC luminosity 
of submm-bump scattered photons. Within the IC picture then, it would be quite natural 
for the variations seen in the L'-band luminosity to be due to ffuctuations in the magnetic 
field. The IC luminosity of submm bump photons, not dependent on B, would remain 
unaffected by such fiuctuations. If relativistic effects are also taken into account (Doppler 
boosting in particular) then the magnetic ffeld ffuctuations we are talking about are actually 
the ffuctuations in the magnetic ffeld of the observable region at any given time. 

However, if the magnetic ffeld were as high as that found in the SED modelling section 
for the IC model {B ~ 210 G), then there would be a cooling break below IR frequencies 
and this picture could no longer work, since above the cooling break the lightcurve traces 
the pure rate of energy injection (no longer with any dependence on magnetic ffeld). We are 
faced with a dilemma (in addition to the small size of the quiescent region which is required) , 
where high magnetic ffelds are required to ffnd an acceptable solution for the observed SED of 
the ffare, but low magnetic ffelds are needed to keep the cooling break above IR frequencies, 
and we can not have both at once. 

If we now examine the SSC scenario, the SSC luminosity goes quadratically in the 
quantity NO"^, which means that a synchrotron lightcurve has a natural width that is a factor 
-\/2 times the width of its SSC emission (for example, if the lightcurves can be described by 
Gaussian proffles f{t) oc exp(— /c(t — to)^/u^^)- Thus a longer duration synchrotron lightcurve 
is expected in the synchrotron case. 

The observations of substructure however are not naturally explained. Within the SSC 
scenario both the SSC luminosity and the synchrotron luminosity are proportional to B^, 
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and fluctuations should thus induce variations of similar strength in both L'-band and X-ray 
lightcurves. As discussed in Section 13.31 substructure would be distinguishable in the X-ray 
lightcurve if it were of similar amplitude to that of the L'-band lightcurve. 

Finally, there is the synchrotron scenario with a cooling break. Understanding simul- 
taneous light curves in this scenario is more sensitive to the time dependent evolution of 
the electron distribution itself and thus self-consistent time-dependent modeling is required. 
We can make at least some qualitative predictions for the light curves of this model: as far 
as the difference in light curve widths is concerned, it seems it would be necessary for the 
cooling break to increase in frequency during the flare which would, as for the submm IC 
case, require a decrease in magnetic field during the flare (see Equation [3]). In contrast to 
the submm IC case, for this case the decrease in magnetic field could occur together with 
plausible values for the magnetic field. 

It also turns out that obtaining substructure in the IR light curve at the same time as 
producing a smooth X-ray light curve could be quite natural in the cooling break synchrotron 
model. Below the cooling break the emitted synchrotron spectrum is sensitive to variations 
in the magnetic field which is expected to be clumpy; while above the cooling break on the 
other hand the synchrotron emission traces rather the rate of energy injection alone. This 
could be expected to be rather smooth. 

We think these aspects add very much to the plausibility of the synchrotron scenario as 
a viable mechanism for the production of the NIR/X-ray flare. 



6. Discussion 

6.1. Why inverse Compton scenarios don't work 

Here we give some analytical arguments that clarify which essential features of our 
multi wavelength observations lead us to exclude the inverse Compton processes as possible 
emission scenarios. 



There are three equations (see lRybicki fc Lightmanlll986l ) which essentially describe all 



important relationships between seed, synchrotron and scattered spectra involved in one 
inverse Compton scattering process. The first of these is the relationship that describes the 
shift in frequency of a seed photon upon encountering an electron of energy 7: 



^IC = I'^^seed 



For the submm IC case, photons are scattered from the submm bump, at Useed ~ 10^^ Hz. 
For the X-ray slope to be soft {ujc < 10^^ Hz) this equation restricts the electron energies 
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of the electrons involved in the flare: 

7 < 1000. 



Secondly, the frequency at which these electrons are themselves emitting synchrotron 
emission is dependent upon 7 and B B 

i/c = 4.2 X lO^B-f"^. (7) 



Since we have already made a restriction on the electron energies, then if we require 
Uc > 10^^ Hz (for a hard uL^, spectral index in the IR), then we find that the magnetic field 
is also restricted: 

5 > 25 G. 

Thirdly, there is the equation relating the ratio of IC to synchrotron luminosity: 

compt ^phjSeed /q\ 

7 " ~Ti — ^ ^ 

-t^synch 

where f/ph,seed is the energy density of seed photons, f/ph,sccd — -f'ph.seed/cA, and A the surface 
area of the electron population producing the seed spectrum. 

In the submm case Lp^^^eed is -^^submm, and the equation can be rewritten for Rq = 
y/A/An as 

With this we obtain an uncomfortably low constraint on the size of (the most luminous 
part of) the quiescent region of Rq < O.lRs, similar to the small sizes we found were required 
in the SED model fitting of Section |H 

We can repeat the above series of arguments for the SSC case to obtain 



7 


< 


100 


B 


> 


2400 G 


Rp 


< 


0.002 R. 



^Note that the formula we give here is for the case where electrons are spiralling exactly perpendicular to 
the magnetic field direction; i.e. the above equation actually contains a sin 9 term {9 the pitch angle of the 
electrons to the magnetic field) which is at maximum 1. If we took a smaller constant value (more realistic 
for an electron distribution with isotropically distributed pitch angles) then the results which follow would 
be even more restrictive. 
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A constraint on the size of the flare emitting region itself also enables a constraint on 
the density of the flare emission region: 



n, ^ 2.4 X 10^ (^yY^^f cm-3 

and we obtain that Ue > 10^° cm^^. As with the density we found in the SED modelling of 
Section HI this is a very high density, several orders of magnitude higher than the density 



inferred for the inner regions of the accretion flow around Sgr A*(~ 10 cm : lYuan et al. 



20031). 



As for the submm IC case, the driver for the high magnetic fields and densities in the 
SSC case is the restriction on the electron energies required by the soft X-ray slope, which in 
combination with a positive IR slope forces the magnetic field to very high values. Adding 
to that the high luminosity ratio, Lx/Ljfj, the size of the flare region is driven to very small 
values which in turn forces the density to very high values. 

These arguments show that for both cases there are three main properties which in 
combination are driving (i) in the submm IC case, the magnetic field to higher values than 
i? ~ 30 G and the IC region to small sizes, and (ii) in the SSC case, to extremely high 
magnetic fields and electron densities. These are: 

1. the soft X-ray spectral index 

2. the hard MIR-IR spectral index 

3. the high luminosity ratio, Lx/Lir. 



6.2. Comparison with past mult i wavelength studies 



Past multiwavelength observations of Sgr A* have favored models where the X-ray 
emission is due to inverse Compton scattering processes, with various possible combinations 
of seed photons and seed electrons from those producing the quiescent (submm) and flaring 
(IR) states. The obvious question is: why do we not find the same? 



Eckart et al.l (120041 ) modelled both IR and X-ray flares as SSC emission. However, that 



the flares at IR wavelengths are SSC emission is ruled out since the polarization of the IR 
flares points to a synchrotron origin. 



.Ecka.rt et al.l (l2006d ) modelled the X-ray flare as SSC emission via the prescription of 
Marschei and in their model the X-ray emission had a hard vLy index of /3 = 0.4 to 
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match their X-ray observations. A hard vLy index in the X-ray will provide no constraint 
on the energies of the electrons producing the IR flare, which explains why these studies 
found this model to be acceptable with reasonable sizes for the quiescent state of Sgr A*. 
The X-ray flares of the ir study were howeve r quite weak and the photon index may not have 
been well constrained (jPorquet et al.ll2008l . Mascetti et. al. 2008 (submitted)). 



A similar story appli es to other invest igatio ns of the X-ra y flare s with IC and SSC 
processes. For example, in iLiu et al.l (j2006bl ) and iMarrone et al.l ( l2008l ) the X-ray emission 
was modelled with SSC, bu t again this was for a hard spectral index in pL^. For these cases 
as for lEckart et al.l (l2006d ) there was no restriction on the electron energies since the X-ray 
flare was taken to be hard, and accordingly it was possible to find a model with reasonable 
physical parameters. 



Yusef-Zadeh et al.l (120061 ) suggested two scenarios where the X-ray emission could be due 
either to (i) submm photons upscattered by electrons producing IR emission, the scenario we 
considered in Section 14.11 or (ii) IR (flare) photons upscattered by the electrons producing 
the quiescent state of Sgr A*. 



The first case is the case we explored in Section 14. 1[ In lYusef-Zadeh et al.l (120061 ) 
however, the treatment involves a power law distribution of electrons and the differential 
flux is calculated for the corresponding power law section of the inverse Compton scattered 
spectrum only. Since the X-ray spectral index in this model was considered to be hard in 
vLy {(3 = 0.4), we have the same situation as with the previous cases; there is no restriction 
on the electron energies producing the IR flare. Therefore, with electrons up to 7 ~ 6000, it 
was possible to find a model that worked using a reasonable quiescent region size {Rq ~ 10 

Rs). 

The second case is a scenario we did not consider in our modelling. In fact it can 
be shown that the inverse Compton luminosity via this process (IR photons scattered by 
submm-emitting electrons) can never exceed the luminosity of the IC case we considered 
previously (submm photons scattered by IR photons). 

For the case of IR photons scattered by submm-emitting electrons, the luminosity is 
(here using u to denote a photon) 



LjciJR i^, submm e" 



2RfL 



submm^ I R 



RlcB^ 



where we have made use of the fact that only some proportion of the quiescent electrons 
(those within the flaring region where the photon density of IR photons is highest) are 
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available to inverse Compton scatter photons: 

T - 

-'-'synch j-,3 ^ submm- 

For the submm IC case that we already considered we have 



Licisuhmm i/, IR e ^ - „ ^9 

KqCB^ 



and the ratio of the two is 



Ljc{lR ly, submm e ) Rp 
L/(7 (submm u, IR e~) Rq 

Thus the X-ray luminosity provided by the IR seed photon case is always going to be, 
at best, comparable to the X-ray luminosity produced in the submm seed photon case and 
will never dominate the emission. In Section 14.11 we found it was difficult to find a solution 
that did not involve an unrealistically small size for the quiescent region of Sgr A*. The 
contribution of IC emission through the scattering of IR seed photons by submm-emitting 
electrons is at best comparable to this emission and can not satisfy the observations either. 

We also compare expectations for the inverse Compton scattered vs. synchrotron 
lightcurves for this scenario. The different widths of the lightcurves is as difficult to un- 
derstand as in the case of inverse Compton scattered submm bump photons; again, the only 
way out may be for the magnetic field to decrease during the fiare, and to be restored at the 
end of the flare. Additionally, if the IR flare provides the seed photons for the X-ray flare, 
then we should expect to see substructures in the X-ray lightcurve of the same order as and 
simultaneous with those in the IR lightcurve. Fluctuations in the magnetic fleld could not 
help this scenario because the IR flare is directly providing the seed photons for the X-ray 
flare in this case. Therefore we conclude that the time-resolved features of the lightcurves 
do not support this emission scenario either. 

Overall, we have covered all reasonable conceivable inverse Compton scenarios for the 
origin of the X-ray flare simultaneous with our IR flare and have concluded that none of 
these inverse Compton scattering scenarios are viable. 

Of past investigations, the most similar synchrotron mo del to those we pre sent in this 



paper is the one-component synchrotron model presented by lYuan et al.l (12004J . see Figure 
3). To our knowledge, this is the only previous work to suggest that both IR and X-ray 
flares are produced by synchrotron emission. Interestingly, in that paper, a one-component 
model for the electron distribution is ruled out on the basis that the X-ray spectrum is too 
soft (r ~ 2) to be consistent with the X-ray flares observed by Chandra. 
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6.3. The electron energy distribution 

In Section |4] where we investigated inverse Compton scenarios for the simultaneous 
IR/X-ray flare (icmodel and sscmodel) we assumed a thermal distribution for the energy 
distribution of the relativistic electrons. The reason for this choice was for ease of calculation 
since the synchrotron emission of a thermal distribution of electrons is described by well 
known formulas (see Appendix). The thermal distribution is also an expected result of 
turbulent heating and ra diative cooling pr ocesses and has been used elsewhere in models 



for the flares of Sgr A* (iLiu et al.l l2006bl ). Nevertheless, one might worry that through 
our assumption of a thermal distribution of electron energies we have limited ourselves to 
a special case and that there may exist other distributions of electron energies that allow 
the inverse Compton scenario to be an explanation for the April 4 IR/X-ray flare after all. 
Although we do not expect our results to be particular to the thermal distribution given that 
the arguments we developed in Section 16.11 were not specific to any particular arrangement 
of the electrons, we would like to here corroborate this expectation through modeling. 

To investigate the senstivity of our results to the form of the electron distribution we 
investigated power law models for the NIR synchrotron emission with the inclusion of inverse 
Compton scattering. We take power law models of the form ^(7) oc 7"^ between energies of 
7mm and ■jmax and calculate the synchrotron and inverse Compton emission by integrating 
the synchrotron spectrum of a single electron of energy 7 over all electrons in the distribution. 
We fixed the parameter 7mm = 20, assuming that the electrons are accelerated out of the 
population of electrons (7 ^ 10 — 20) that create the quiescent state synchrotron emission, 
but we also found that the results are not at all sensitive to the value of 7mm as we show 
in Figure [H] where we show also a model with 7mm = 1- As we did for our previous set of 
models, we fit the calculated SEDs to the observed data using XSPEC, naming this model 
powerlawicmodel. 

We tested three models, with p = 2, p = and p = —2, corresponding to falling, fiat 
and rising electron distributions respectively. In Figure E] we show the emitted SED and 
the electron distributions that correspond to the best fit in each case. The distribution of 
electron energies is irrelevant to finding a good fit to the observed data, and the SEDs for the 
best fits for all models appear similar. Amongst the best fit models the electron distributions 
show a lot of variety. For a power law of electrons with index p = —2 to which we also add 
an exponential cutoff, we find that the electron distribution of the best fit approaches the 
thermal distribution corresponding to the best fit for icmodel, as we would expect. 

Interestingly, for this variety of models, the electron distributions all meet at a particular 
energy, 7 ^ 1000. It appears that this energy is important {even, the only important energy) 
for the submm IC scenario, and we note that it corresponds exactly to the upper limit on 
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Fig. 8. — Upper panel: Best fit SEDs where the NIR synchrotron emission is produced by a 
power law distribution of electrons with index p = 2, 0, —2 (magenta, green and red dashed 
lines). The best fit model for a thermal distribution of electrons is shown in solid blue. 
All models appear relatively similar. Note that the (model-dependent) X-ray data is only 
shown for the thermal distribution of electrons case, but it should be similar for the other 
models. Lower panel: The electron energy distributions of best fit submm IC models each 
assuming a different underlying form for the energy distribution. We compare power law 
models (n(7) oc 7~^) of index p = 2, and -2 (magenta, green and red respectively). We 
also demonstrate that in a power law model of index -2 with an exponential cutoff (blue) 
we found the same energy distribution as in the thermal distribution case (in black; note 
that we used analytical equations to calculate the synchrotron emission for the thermal case 
in icmodel and sscmodel of Section H]). Although the electron distributions all look very 
different, they all cross at an energy 7 1000. 
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electron energy that we found in Section 16. ![ for the submm IC case. 

In Table |2] we give parameter values for the power law model with p = 2. We find 
that in general the parameters, in particular the confidence intervals, are similar to those in 
Table [1] for icmodel. Of particular interest is to see whether a larger value of Rq can be 
accommodated to revalidate the inverse Compton scenario. The upper limit (90% confidence) 
for this parameter however of -Rq < 0.32 however, is still very small and does not allow for 
the this scenario to be compatible with the size measurements of Sgr A*. It appears that our 
finding for the submm IC scenario, where we find the implied size for the submm-emitting 
region must be much smaller than the observed size, is robust and is not sensitive to the 
form of the electron distribution. 

Beyond assuming a particular form for the electron distribution, we should keep in mind 
that in reality the electron energy distribution is an evolving function that depends on the 
details of injection/escape and electron cooling. A truly self-consistent approach would allow 
the electron distribution function to be determined from the parameters themselves (such 
as, e.g. the magnetic field, which determines how long the electrons take to cool and thus 
has an important effect on the shape of the electron spectrum). The synchrotron spectrum 
would then be calculated directly from the time-dependent energy distribution function. 

As an example of why the electron distribution should be self-consistently calculated 
from the model parameters, we note that both of our inverse Compton models require high 
magnetic fields. Even for the submm IC model with the lower of the two magnetic fields, 
the cooling timescale (Equation [1]) is shorter than the dynamical timescale already at IR 
frequencies. The electron distribution in the case of such high magnetic fields would then 
be already dominated by cooling at L'-band wavelengths. In essence, a further requirement 
on our models is that we hope to find a solution with B < 30 G in order to keep the 
cooling break above the NIR band. However, as we have found in our analysis, there are 
certainly no solutions for the SSC model with such low magnetic field, and the submm IC 
only accommodates such a low magnetic field at the edge of its 90% confidence region. 



6.4. The jet model 



The model spec trum we used to calculate the photon density of submm photons in 
icmodel was that of lYuan et al.l (120031 ). However, this is not the only model that describes 
the qui escent spectrum of Sgr A*; the quiescent spectrum is also well described by a jet 
model JPalcke fc MarkoflJboOoh . 



To test whether our conclusions were robust with respect to the assumption of a qui- 
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Table 2. Models: Fit Parameters 



Fit to Mean Fluxes & NIR, spectral index 
Parameter powerlawicmodel 



Nh [x1022 cm-2] 


11.8 


(10.0, 14.3) 


B [Gauss] 


330 


(22,2100) 


^max [ ] 


1100 


(280, 1600) 


A^e [xlO'"^ electrons] 


20 


(1.7, 1900) 


Rq [Rs] 


0.028 


(0.00003,0.32) 


Rf [Rs] 


> 0.02 




P 


2.0 




X2 / d.o.f. 


70.3 /75 




reduced 


0.94 




Violates 3cr MIR upper limit? 


No 





Note. — Summary of best fit parameters for 

powerlawicmodel fixing p = 2. The constraints on the 
parameters are very similar to those found for the icmodel of 
Section |4] and appear to obey the constraints of Section 16.11 




Fig. 9. — A submm inverse Comp ton model for the X-ray flare using a different model (the 
jet model, iFalcke fc Markoflll2000h as template for the submm photon spectrum. The dashed 
gray line shows the jet model which is used as the source of seed photons. The solid blue 
line shows the best fit XSPEC model. 
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escent model spectrum we created a new XSPEC m odel in which the submm spectrum to 
be upscattered was that of iFalcke fc Markofj (120001 ) . Figure [9] shows the best fit model for 
this case, and the best fit parameters are very similar to those of icmodel, with somewhat 
higher magnetic field {B ^ 300 and lower electron temperature 6e ~ 80). The size of the 
quiescent region, Rq is still very small {Rq ~ 0.02) for the best fit model. Thus we can 
conclude that our findings are also not sensitive to the particular model we used to model 
the submm photon density. 



6.5. Substructure in the context of an orbiting hot spot model 



It is worth making a closing comment on our results in the context of the orbiting blo b 
model JCenzel et al.lboosi : ISroderick fc Loeblboosi : iMever et al.]l2006al : iTrippe et allboOTh . 
where the substructures observed commonly in the IR/NIR lightcurves are postulated to be 
due to relativistic beaming on the approaching side of the hot spot's orbit about the SMBH. 

Due to the fact that the X-ray lightcurve is so smooth, and the substructure only shows 
up in the NIR lightcurve, then naively, the hot spot model does not seem to be compatible 
with the observed lightcurves; if relativistic beaming is occurring, it is not obvious how such 
structures could not also show up in the X-ray lightcurve. 

We think this aspect of the multiwavelength emission from Sgr A* certainly deserves 
consideration in models describing the flare emission from Sgr A* via a hot spot model. The 
different widths of the lightcurves are an additional important clue in this regard. Conversely, 
it is also obvious that given its likely proximity to the SMBH of ~ 4 x 10^ Mq relativistic 
effects can have a considerable influence on the observed emission from Sgr A* and should 
be taken into account in models that try to explain multiwavelength properties of flares from 
Sgr A*. 



7. Conclusions 

We have presented the results of a simultaneous multiwavelength campaign at L'-band, 
X-ray and MIR wavelengths carried out on April 4 2008. We summarize the main observa- 
tional results as follows: 

• The L'-band and X-ray flares were simultaneous to within 3 minutes. 

• The L'-band flare is much broader overall than the X-ray lightcurve. 
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The L'-band flare showed significant substructure with a timescale of ~20 minutes, 
while the X-ray fiare showed no corresponding significant substructure. 

The vLy spectrum increases between 11.88 /xm and 3.8 /im. 



The X-ray fiare was very bright and soft in z/L,^ {(3 = —0.3 ±0.3; iPorquet et al.ll2008l ). 
The emission region must be small, < 1.5Rs- 



We have drawn conclusions about the emission mechanism behind the fiares which can 
be traced essentially to the hard MIR-IR uLi, spectral index, the soft X-ray z/Lj, spectral 
index, and the high ratio of X-ray to IR luminosity. Our conclusions are: 



• We strongly disfavor the SSC case due to the extremely high magnetic fields and elec- 
tron densities required to reproduce the observed data. Both quantities are several 
orders of magnitude larger than the values expected for the inner regions of the accre- 
tion fiow about Sgr A*. 

• We disfavor the submm IC case due to the high submm photon densities required to 
produce the observed X-ray emission, which imply a quiescent region size much smaller 
than implied by the size measurements. The IC case is also only marginally compatible 
with physically plausible magnetic fields of 10 — 30 G. 

• We also disfavor a different IC case (where the X-ray flare is due to IR seed photons 
upscattered by quiescent state submm-emitting electrons) since the inverse Compton 
scattered luminosity in this case is always dominated by the inverse Compton scattered 
luminosity of submm photons scattered by IR-emitting electrons. 

• We favor a synchrotron scenario where the emitted spectrum flattens towards X-ray 
energies due to a cooling break. This scenario can be achieved with physically plausible 
magnetic fleld strengths and holds promise to explain more detailed structures of the 
light curves. 



• For both the synchrotron and IC cases, the relatively shorter duration of the X-ray 
flare and perhaps the substructure of the IR emission is most plausibly understood as 
a result of a transient decrease of the magnetic fleld in the region of the flare. Such a 
decrease may point toward flares being triggered by conversion into electron heating 
of stored magnetic energy, such as a magnetic reconnection event. 
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A. icmodel and sscmodel 

Both icmodel and sscmodel are subsets of a single model with five parameters B, N, 
Rf and Rq. 

For this model, we assume a spherical homogeneous emission region (containing the 
transiently heated/accelerated electron (fiare) population emitting in the IR) of size Rp and 
of electron density Ue = N/{47iR^). The region contains a homogeneous magnetic field 
of strength B. We choose to model the electron distribution with the thermal electron 
distribution (i.e., a relativistic Maxwellian) 



n(7) 



ne7Vl-l/7^ 



with 6e = kTe/rrieC^ the dimensionless electron temperature, and K2{x) a modified Bessel 
function of the second kind. The thermal distribution is a good example of an electron 
distribution of a c haracteristic ener gy (i.e. 7 ~ 9e), and is an expected result of turbulent 
heating processes (ILiu et al.ll2006bl ). 



The emission coefficient is approximated by (IMahadevan et al.lll996h 



Ju,th 



2u 



This preprint was prepared with the AAS I^TJt;X macros v5.2. 
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in units of erg s ^cm ^Hz ^ster ^, and with 

M(x) = 4.0505ax-^/^ exp(-1.8896x^/^)(l + OAObx'^/^ + 0.5316cx^^^^). 
The absorption coefficient is 

Oiu,th = ju,th/B^{Te) = j„,th -^j^ (exp (hu/kTe) - 1) 

and the resuhant synchrotron spectrum is computed for our spherical flare region (including 
optical depth effects) by the equation of radiative transfer 

Lj,s = 47r / — ( 1 — exp{—a,y^/Rf^'^-r^] Airrdr 
Jo ^ ' 

which for small optical depth (an optically thin flare) simplifies to 

'47r 



4vr —R% Ju. 



We c ompute the Inverse Compton scattered luminosity in our model through (IBlumenthal fc Gould 



197(1 1 

where 7, e and ei = hu are the electron energy, initial and scattered photon energies respec- 
tively, and the quantity 

dN^Jdtdeide = Sarc nph{e)/A^^e [2glng + (1 + 2g)(l - g) + 

where = Ae'j/mc^ is the Compton factor, q = ei/Veipfmc^ ~ ^i)- The photon density 
nph{e), of flare state photons is determined from the model flare's luminosity by nph{e) = 
nph{^)/h = Ly/i^ATih^vcR^). Instead of modell ing the qui e scent state ourselves, we use 
the luminosity given by the quiescent model of lYuan et al.l (120031 ) . which reproduces the 
observations quite well. Using this quiescent model as input spectrum L;^ q then, the photon 
density of quiescent state photons is nph{e) = L„^q/ {A/STch'^ucRQ'^). 

We implement icmodel and sscmodel via the above prescription, with only one differ- 
ence concerning the parameters Rq and Rp- For icmodel, Rp is not an input parameter, 
and is merely set to Rp = Rq (this simply ensures that a minimum of SSC emission is pro- 
duced, so that a solution where the X-ray emission is due to inverse Compton scattering of 
submm photons may be found). Similarly for sscmodel, Rq is not an input parameter, but 
is set to an arbitrarily high value (in this case to ensure that very little inverse- Comptonized 
submm emission is produced and a solution is found where the X-ray emission is due to SSC 
of IR/NIR photons). 
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B. powerlawcool 

This model can be written as 

uL,, oc < 



Z/(2-P)/2 Z/c < 7 < I/n 



The coohng break, Vf. (Equation [3]) occurs at = 2.56 (5/30 G) ^ x 10^^ Hz, and the 
index p corresponds to the particle index in the underlying electron distribution: 



n(7) oc 



7 ^ 7min < 7 < 7c, 

7"^^^^^ 7c < 7 < 7max. 



